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Abstract 

We perform a complete calculation of the relic abundance of the KK-photon 
^) ', LKP in the universal extra dimension model including all coannihilation channels 

and all resonances. We show that the production of level 2 particles which decay 
^ ' dominantly into SM particles contribute significantly to coannihilation processes 

Q , involving level 1 KK-leptons. As a result the preferred dark matter scale is increased 

to = 1.3 TeV. A dark matter candidate at or below the TeV scale can only 
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be found in the non-minimal model by reducing the mass splittings between the 
KK-particles and the LKP. The LKP nucleon scattering cross section is typically 

10 



(— I I small, (T < 10 pb, unless the KK-quarks are nearly degenerate with the LKP. 

& 1 Introduction 

One of the most attractive explanations to the dark matter(DM) problem is a new weakly 
, interacting particle (WIMP) present in extensions of the standard model. Supersymmetry 
t — I and extra dimension models are the leading candidates for physics beyond the standard 
^ ■ model that also propose a WIMP dark matter candidate. Among the extra dimension 
^ ■ models, the UED scenario [1] where all standard model particles are allowed to propagate 
^ . freely in the bulk is of particular interest. In this model momentum conservation in the 
extra dimensions entails conservation of KK number. Orbifolding is required to obtain 
chiral zero modes from bulk fermions, and breaks extra dimensional momentum conser- 
vation. However, there remains a discrete subgroup, KK parity, thus the lightest KK-odd 
particle is stable. In the minimal universal extra dimension model (MUED) the dark mat- 
ter candidate is in general a vector particle, B^, the Kaluza-Klein (KK) level 1 partner 
5^ i of the U(l) gauge boson. One characteristic feature of the model is that the dominant 
annihilation channels of B^ are into lepton final states. The observations of PAMELA [2] 
and Fermi [3] hinting at an excess in the leptonic channel with no counterpart in the 
antiproton channel has therefore renewed interest in the UED models. Note however that 
this excess can be explained with astrophysical processes and does not necessarily require 
a dark matter interpretation [1]. 

In the MUED model all KK states of a given level have nearly the same mass at tree- 
level, n/R, where R is the size of the compact dimension. The mass degeneracy is lifted 
only by standard model (SM) masses. Radiative corrections at the one-loop level further 
induced mass splittings among level 1 particles, these mass splittings are however small 
for all weakly interacting particles. This means that co-annihilation channels naturally 
play an important role in the computation of the relic abundance of dark matter. The 
relic abundance of B^ in the MUED model was first computed in [5] and included the 
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coannihilation of SU(2) singlet KK leptons, the Next-to-Lightest KK particle (NLKP). 
Apart from the mass difference entering the Boltzmann equation, this calculation was 
performed in the limit of degenerate masses for all particles at a given KK level. It was 
shown that the effect of KK-leptons is to increase the dark matter relic abundance. This is 
because the coannihilation channels are not as efficient as the main annihilation channels, 
the increase in the number of effective degrees of freedom then induces an increase in 
the relic abundance contrary to what usually occurs in the MSSM. When the light Higgs 
mass is large (around 200GeV), the KK-Higgs is the NLKP, coannihilation processes with 
first KK level Higgs particles lead to a decrease of the relic abundance [6] increasing the 
preferred mass scale for dark matter. The relic abundance including all coannihilation 
channels as well as a precise evaluation of the KK masses were later computed in [3 [S] . 
The impact of the precise value of the mass splittings between the particles of the first KK 
level and the LKP was also analysed by going beyond the MUED framework and treating 
the mass splittings as free parameters. Note that in non minimal UED versions, it is 
possible to modify the mass splittings by making different assumptions on the boundary 
terms at the cut-off scale. 

At tree-level the second level KK particles have masses twice as large as the ones of the 
first KK states. It is therefore natural to expect an enhancement of annihilation channels 
of 1st level KK particles due to the exchange of a s-channel 2nd level KK particle near 
resonance. For B^B^ annihilations, the only potential resonance is the KK partner of the 
Higgs, h"^. This particle either decays into B^B^ or into standard model particles. It was 
shown in [S] that the loop-induced decays into standard particles dominate. The inclusion 
of h"^ in the annihilation processes thus reduces significantly the relic abundance O ITO] . 
Other resonance effects occur in coannihilation channels, for example e^e^ — )■ Z'^. The 
computation of the relic abundance including all coannihilation channels and level 2 reso- 
nances from the Higgs and gauge boson sector was performed in [TT]. The effect of level 
2 fermions was not included as those particles do not decay dominantly into SM particles. 
There is however an additional effect from level 2 particles that has been overlooked until 
now, it is the production of a KK-even level 2 particle in association with a SM particle 
in the final state. This level 2 particle then decays through loop induced processes into 
standard model particles. For coannihilation processes where channels with B^ in the 
final state are usually kinematically accessible, there can be strong enhancement due to 
resonance effects from a variety of level 2 particles, for example the coannihilation process 
e^i?^ — 7- — eB^ can be strongly enhanced by the contribution of the level 2 lepton. The 
coannihilation processes can even dominate over the annihilation channels thus reducing 
the relic density oi B^ . 

In this paper we compute the relic abundance of the LKP including all coannihilation 
channels, all possible level 2 resonances as well as level 2 particles in the final state. We 
take into account electroweak symmetry breaking effects that were neglected in previous 
calculations. These will for example impact the masses of the KK-Higgs particles. We 
consider the case where the DM candidate is the KK partner of the photon, 7^, rather 
than B^. Note however that the LKP is dominantly compose of B^ with only a small 
component. The new contributions from annihilation channels into gauge bosons 
will be suppressed due to this small mixing angle. Finally we include couplings between 
level-2 gauge or Higgs bosons and standard model fermions. In particular we perform a 
one-loop computation of the couplings of the level 2 KK-Higgs into quark pairs updating 
previous calculations [TT] by including the contributions from n = 1 scalars and weak 
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gauge bosons. We work first in the context of a minimal UED modefl then in an extension 
of the minimal model by treating the mass corrections in the fermion sector as arbitrary 
parameters. This is done to illustrate the strong dependence of the rehc density on mass 
corrections to the KK states. Indeed the mass corrections influence the coannihilation 
suppression factor, determine how strongly the channels with level 2 KK-particles in the 
final state are kinematically suppressed as well as whether the (co)-annihilation process 
benefits from a resonance enhancement. In the MUED model we find that the mass 
splitting is such that the second KK level in the final state play an important role. Thus 
the relic abundance is reduced and the preferred value for the DM mass, consistent with 
cosmological measurements, Qh"^ = 0.1120 ± 0.0056 [13], is shifted above the TeV scale. 
We also give predictions for direct detection rates and for KK particles production at the 
LHC. 

The paper is organized as follows: in section 2 we briefly present the UED model and 
give explicit expressions for mass splitting of level 2 particles and loop induced vertices 
between level 2 and SM particles. In section 3 we discuss DM (co-)annihilation channels 
and DM observables. In section 4 the results for the relic abundance and the detection 
rates are discussed. Section 6 contains our conclusions. 

2 The UED model 

In the UED model all SM fields are upgraded to higher dimensional fields that propagate in 
the bulk of the flat and finite spatial extra dimensions. Among many space-time structures 
that can reproduce the SM as a low energy effective theory, we focus on the simplest 
construction: the five- dimensional space-time with the extra dimension compactified on 
S1/Z2. Integrating out the extra dimension leads for every 5D field to an infinite tower 
of 4D Kaluza-Klein modes. In the limit of 5D Lorentz invariance, the mass spectrum of 
the n^^ KK modes is n/R, where R is the radius of the compactified dimension. The Si 
compactification destroys the 5D Lorentz invariance, and the Z2 orbifolding violates 5D 
momentum conservation. The dispersion relations of 5D particles are modified, resulting 
in shifts of the KK mass spectrum from n/R and in new decay patterns. The Z2 orbifolding 
leads to chiral zero modes (the SM fermions) and removal of the zero modes of the extra 
dimensional components of the gauge bosons, and allows for 4D brane interactions on the 
orbifold fixed points. Even though the Z2 orbifolding violates momentum conservation in 
the extra dimension, we can retain a parity symmetry under the flip of the fifth coordinate. 
From the four-dimensional viewpoint, this Z2 parity prevents mixing between even and 
odd KK levels, thus guaranteeing the stability of the lightest KK-odd particle, the LKP. 
The coupling of the KK modes are generated from the 5D Lagrangian after expansion 
over the KK modes. 

We first construct a gauge invariant 5D SM Lagrangian consistent with the above 
observations. Kinetic terms violating the 5D Lorentz symmetry are compatible with the 
5D gauge invariance, and can practically produce the loop-induced mass shifts of the 
MUED model, as we will see shortly. We consistently implement electroweak symmetry 
breaking: the vacuum expectation value of the Higgs scalar field is involved not only in 

^Note that our version of the minimal UED model differs slightly from the MUED model of Ref. [T^ . 
in particular as concerns the masses of the n — 2 gauge bosons and the couplings of level 2 Higgses to 
SM fermions as detailed in section 2. 
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the KK electroweak gauge boson and KK fermion mass matrices, but also in the KK 
Higgs boson mass sphtting. In addition, mixing angles caused by electroweak symmetry 
breaking are taken into account in any interaction term. Gauge-fixing functions and 
ghost terms are introduced using 5D Goldstone and ghost fields. As far as the bulk 
5D Lagrangian is concerned, the KK number is conserved at each vertex. Then, for 
phenomenological studies, we adjust the KK mass shifts to the MUED ones, which are 
radiatively generated, and add to our UED model KK-number-violating direct couplings 
between level 2 particles and SM particles as perturbation because such couplings are 
necessarily induced at the loop level. The complete description of the Lagrangian and of 
the gauge fixing procedure is provided in another publication [H]. Here we only describe 
the mass spectrum including radiative corrections as well as loop induced decays of level 
2 particles into pairs of SM particles. 



2.1 Mass corrections 

The UED model is an effective field theory valid up to a cutoff scale A. At loop level mass 
shifts are introduced by bulk corrections and by localized brane terms. Bulk corrections 
affect only the gauge bosons KK states while localized brane terms induce mass shifts 
for all particles. In the minimal UED model one chooses vanishing boundary terms at 
the cutoff scale. Nevertheless after running down to the electroweak scale radiatively- 
generated bulk and brane terms affect the mass spectrum lifting the degeneracy of the 
KK levels. The mass corrections to 5" and were computed in Ref. [12]. The squared 
mass matrix for neutral gauge bosons reads 



(n/ R)"^ + m\c\^ + 5m^n —m\cwSw 



where sw-, cw are the sine and cosine of the weak mixing angle and 

e 2 5 5'C(3) 1 Ihg'n^ 1 , 

^"^^^ = "2 16^^ + ^:R^ 16^ [j- ) 

The renormalization scale /i is chosen to be the mass scale of the 1^* KK mode, /i = -R~^, 
and g,g' are the SU(2) and U(l) gauge couplings respectively. The mass of the LKP, 7^, 
is obtained after diagonalisation of the 5^, mass matrix. 

Zf) = W^,3Wcos<)-i?(-)sin<\ (4) 

where is the n}^ level mixing angle which is near for n > 1. The mass of the LKP is 
given approximately by the mass of in Eq. [T]and deviates from at the permil level 
or less. To take into account the loop corrections to the mass terms in a gauge invariant 
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manner in the 5D Lagrangian we introduce a correction to the fifth components of the 
SU(2) and U(l) gauge fields, namely, 

= -\b,,B^^ + ^B^,B^ + ^W^.W^'^ + ^W^5W5^ (5) 

This improved tree-level Lagrangian reproduces the loop-corrected mass terms once 
the parameters Z^, Z\y are fixed so that the new mass matrix matches Eq. [H see ref. [Tl| 
for a complete description. The gauge boson mass matrix now reads 

^2(n) ^ j Zwin/Rf + mlcl^ -mlcwSw \ 

As long as rrih is light (near 120GeV), the next to lightest KK particles are the right- 
handed KK leptons whose mass corrections are also governed by the U(l) coupling. The 
mass eigenstates are obtained after diagonalisation of the KK lepton mass matrix, 

j^{n) ^ I n/R + 6mi2 mi \ = 1 ^^l^^^ | (7) 

I m -n/R- 6min j \ rrii -Zi^n/R j 

where the parameters , Zi^ describe the loop- improved masses in the tree- level gauge 
invariant 5D Lagrangian, explicitly 

= ^l^I^D^^Pl - Z/,?i7'^5^L + ^Ril^D^i,R - ZiJ^^^^'D.^R (8) 

where D^^^ are covariant derivatives. In the MUED model, the corrections for the KK 
partners of singlet leptons read 

^"^'"« = 7T7JZ?7^1og ■ (9) 



These corrections give a mass splitting at the percent level between the LKP and the 
NLKP. The mass splittings for the KK lepton doublets and the KK gauge bosons are 
governed by the SU(2) coupling and are at the few percent level while those for coloured 
particles are governed by the SU(3) coupling. For the top quark an additional mass 
splitting is induced by a term proportionnal to the Yukawa coupling. The mass differ- 
ence between the coloured particles and the LKP are generally large enough (15-20%) 
that the Boltzmann factor suppresses the contribution of the coloured states in coanni- 
hilation channels. To implement the mass corrections for KK-quarks we introduce three 
new parameters in the Lagrangian Zq^^, Z^^^ Z^^ and follow the same prescription as for 
leptons. 

In the Higgs sector, the KK charged Higgs /i^" is the lightest particle, the CP-odd 
Higgs a" the next to lightest, and the CP-even Higgs the heaviest at each KK level, 
the masses read 



(10) 



5 



where 

^™5,.^(|/4«'^-A.)^^log(^) (11) 

where the Higgs quartic couphng, Xh is related to the mass m/j = \/XhV. Since the quartic 
couphng induces a negative mass correction to the Higgses, for a large value of rrih (and 
Xh) the charged Higgs can become the NLKP or even the LKP [B]. In the 5D Lagrangian, 
corrections to the KK masses of the Higgs doublet (p are taken into account by introducing 
the parameter 

C^^ = _ Z^D,(P^D,<P - /iVV (12) 

The complete mass corrections in the MUED model are given in Ref. [12]. Note that 
when computing the spectrum we fix the gauge couplings at the electroweak scale. 



2.2 Decays of level 2 particles 

At tree-level the n = 2 KK-particles can decay either into two n = 1 KK-particles or 
into another n = 2 KK particle and a SM one. Both these channels are kinematically 
suppressed so that the dominant decay mode can be a loop-induced two-body process 
into SM particles. In particular the n = 2 partner of the LKP will decay dominantly into 
SM fermions. The coupling of 7^ to standard model fermions is a loop-induced process 
involving triangle and self-energy diagrams with level 1 fermions and a level 1 gauge boson. 
The vertices B'^ff and W^ff have been computed in [12] and are related to the mass 
corrections from the boundary terms. The dominant loop-induced coupling comes from 
the level 1 gluon and level 1 quark exchange diagram in the B^qq vertex. In particular 
the vertex B^ti receives corrections that are proportionnal to the top Yukawa, yt. The 
loop-corrected vertex reads 

+ y„(fg"-+6gi-3vmi+'i!>))tBi (13) 

where Yt^, Yt^^ are the hypercharge of ti and t/j. Note that there is a partial cancellation 
between the gluon and Yukawa contribution. For the light quarks one gets the same 
expression neglecting the Yukawa term while the couplings to leptons are suppressed as 
they receive no contribution from gluon exchange nor from large Yukawas. The width of 
7^ is around 1 GeV when R^^ = 1 TeV and the branching ratios do not vary much with 
R or AR, with 13.9% into ti , 29.7% into cc and uu, and 8.6% into each of the d-type 
quarks. 

The Higgs bosons, h"^, and h"^^, are the other level 2 particles where the dominant 
decay mode is into SM particles. The level 2 CP-even Higgs, h"^, predominantly decays 
into ti through the radiatively generated vertex h'^ti. The t^-g^ one-loop contribution to 
the vertex h'^tt is found in Ref. [Tl]. In addition, we include the vertex corrections from 
the t^-W^ and t^-B^ one-loop diagrams, as well as the contributions from the kinetic and 
mass mixings between the level 2 particles and the SM particles on the external legs. 
The latter contributions stem not only from the gauge interactions but also from the 
top- Yukawa interaction and Higgs self-interaction, which allow the level 1 Higgs bosons 
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to run in the loop. The full expression for the vertex h?tt at the one-loop level is given by 

^ - it^ {''^- - T^' ^ k - ^ ''^) if) '''' ■ (^^) 

The coefficients of the vertices, a'^tt and h'^~tb, have the same form due to gauge invariance. 
In our analyses, we omit the loop- induced h^hW-tjpe vertices as such interactions are 
suppressed by the weak coupling constant. Loop diagrams that involve the Higgs vacuum 
expectation value are also neglected. Finally the loop induced vertices to lighter quarks, 
including bb are negligible. Given the above interactions, the total decay widths of h"^ 
(a2,/i2±) are 395(400,401) MeV for R-^ = 1.3 TeV, nih = 120 GeV and AR = 20. The 
branching ratio of each level 2 Higgs boson into the third generation SM quarks is more 
than 99%. 

Notice that in the MUED model the kinetic and mass mixing terms, which contribute 
to the direct couplings of level 2 particles with SM particles, are proportional to the 
radiatively-induced mass shifts of the level 2 particles. In the scenarios with arbitrary 
mass splittings we discuss later, we use the same KK-number-violating interactions as 
in the MUED model. As long as the arbitrary mass splittings are introduced at the 
tree level with the cutoff scale fixed as assumed in this paper, only the argument of the 
logarithm in the KK-number-violating couplings is affected. Since we use the leading log 
approximation, to take into account such effects is meaningless. 

3 Relic abundance 

To compute the relic density of dark matter, we have implemented the UED model in 
CalcHEP [15] and micrOMEGAs2.4 [ISl H?]. For this we relied on LanHEP and the new fa- 
cilities to project the 5D Lagrangian into a 4D one [18]. Implementations of the MUED 
model in CalcHEP [TS] and in FeynRules [20] are available and aim primarily at col- 
lider studies. Here, we briefly describe differences between the implementation of our 
phenomenological UED model and that of the MUED. First, we have implemented the 
classical 5D SM Lagrangian with the wave function factors, Zi, which allow for arbitrary 
mass splittings provided that the mass corrections of KK particles are proportional to n. 
The Higgs sector and the electroweak symmetry breaking are fully taken into account: for 
example, the mass splitting among the KK Higgs bosons, Eq.(9), is automatically gener- 
ated. For simplicity, we neglected the SM quark mixing angles and the renormalization 
group running of the coupling constants, both of which are included in Ref. [19]. In the 
realisation of Ref. [20] the level 2 KK-particles are not included as well as the charged 
and CP-odd level 1 KK-Higgses. Notice that there exist KK Goldstone and ghost flelds 
corresponding to the higher level massive gauge bosons, these were not considered in the 
preceding implementations. We have implemented both the unitary and the Feynman 
gauges, and checked that the 5D gauge invariance is retained. The loop induced vertices 
of the level 2 KK-particles with SM particles described in section 12.21 are then added to 
the model flle by hand, so that the decay widths of the level 2 particles are automatically 
computed. The model will be described in [T3]. 

To discuss MUED phenomenology, the wave function factors are adjusted to reproduce 
the loop-induced mass spectrum of the level 1 KK particles: the masses of the higher 
KK particles are given by integer multiples of those of the level 1 KK particles. In the 
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MUED model, this holds true for all corrections except the bulk corrections on the gauge 
bosons, Eqs. [2|3| and those induced by a shift in the renormalization scale. In total, 
our approximation leads to a shift on the mass of level 2 KK gauge particles less than 
per-mil for 771^2 and 1.2% on the mass of W^,Z^ as compared to the MUED case. In 
comparison, changing the cut-off scale from AR = 20 to AR = 50 leads to 1.6% shifts on 
these masses. Note however that even a small shift for the mass of a level 2 KK-particle 
can give rise to a noticeable discrepancy. For example, in most of the parameter space 
of our model, 7^/1^"^ — )■ h'^'^ — )■ tb gives the dominant contribution to the relic abundance 
among loop-induced resonant processes, while h^^ is light enough that such a pole is not 
reached at the LKP decoupling temperature in the MUED model, as we will see. 

The relic abundance calculation when coannihilation channels are present involve gen- 
eralizing to the effective annihilation cross section. To display explicitly the dependence 
on the mass differences we write (av) as 




where {<Tv)ij is the thermally average cross section for annihilation of any pair of n = 1 
KK particles into KK-even particles and Qi are the number of degrees of freedom. The 
Boltzmann suppression factor is Bi = exp'^'^^^'^f where Tj ^ mLKp/25 and Amj = 

m - TTlLKP- 

When computing the relic abundance we sum over all n = 1 KK particles and include 
all possible SM particles in the final state as well as the KK-even n = 2 particles. Indeed 
as mentionned above both 7^ and the level 2 scalar, pseudoscalar and charged Higgs 
(/i^, a^, h'^^) decay preferentially into SM particles. Furthermore their decay rate, typical 
of a weak interaction process, is much faster than the Hubble expansion rate at freeze-out 
temperatures. These processes therefore contribute to the rate of annihilation of dark 
matter particles. Other n = 2 KK particles, notably KK-leptons, can also be produced 
in the final state in coannihilation processes. These particles can decay into other n = 1 
KK states as well as into SM particles. One can show that the contribution of ra = 2 
KK fermions in the final state to the effective annihilation cross section can be taken 
into account by multiplying their production cross section by the branching fraction of 
KK fermions into SM particles. We have modified micrOMEGAs to take this factor into 
account. Note however that the contribution of KK leptons in the final state is only 
at the percent level so this effect is small. The most general case with large branching 
fractions of heavier level 2 particles into a lighter level 2 particle and a SM particle 
would necessitate modification of the Boltzmann equation for the computation of the relic 
abundance. However such cases are not relevant in the UED model under consideration. 

We have also checked the consistency of our numerical results for the main annihilation 
channels with Ref. [7], good agreement was found when we removed the contribution of 
the level 2 particles coupling to SM particles. Furthermore our results qualitatively agree 
with [9] when we ignore level 2 particles in the final states (but include h"^ exchange). 
Some numerical differences are found due to the more accurate computation of the h'^ 
decay width as described in section 12.21 

The relic abundance calculation can be rather slow in micrOMEGAs because of the 
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very large number of coannihilation channels of the UED model. To avoid computing 
unnecessary coannihilation channels in the calculation of the contribution of individual 
channels to the relic abundance we do not include channels for which the Botlzmann 
suppression factors BiBj < 0.01 (this is done by setting B^ps = 0.01 in micrOMEGAs). In 
the computation of the relic abundance on the other hand we set BiBj < 0.00001. Note 
however that setting B^ps = 0.01 is sufficient to compute Qh"^ with an accuracy below 1% 
as we have tested with several input parameters. 

4 Results 

The minimal UED model contains only three free parameters: R~^, A and rrih. We 
vary the mass scale of KK states in the range 0.4-1.8 TeV. Indirect constraints from 
electroweak precision observables [211 122] and rare processes such as 6 -^■ S7 [231 El] set 
the lower bound between > 0.4 — 0.6 TeV while larger values will lead to too much 
dark matter. We vary rrih = 114.4 — 300 GeV to comply with the LEP limit [23] and to 
include the whole region where the LKP is a neutral particle. For the cut-off scale we 
choose the range AR = 20 — 50. The upper range for the cut-off scale is motivated by 
the computation in Ref. [26] of the renormalization group running of the couplings which 
shows that the U(l) coupling blows up at -E = 50 TeV when R~^ = 1 TeV. 

For rrih = 120 GeV, AR = 20, the value of Qh"^ computed with only the main anni- 
hilation channels shows that 7^7^ annihilate efficiently and that a rather heavy mass is 
preferred {m}^ = 800 — 950 GeV), see Fig. [H The dominant channels do not vary much 
with the scale, for R~^ = 1 TeV they are into II (19% each), it, (21%) qq (17% altogether) 
and uu (3%). The contribution from W,Z boson final states is around ~ 1% and is sup- 
pressed due to the very small B^ — mixing. The enhancement of the tt channel as 
compared to other quarks is due to the contribution of h"^. In Fig. [H we also display the 
value of Qh"^ when neglecting the h"^ exchange in s-channel, this is done by removing the 
loop-induced vertex /t2//- The contribution of /i^ induces roughly a 10% decrease in Qh'^ 
as was found previously [9]. 

Including all coannihilation channels while forbidding the production of 7^ and other 
level 2 particles as well as loop-induced couplings between level 2 particles and SM ones 
leads to an increase in the relic density [S]. This is because the coannihilation cross 
sections (predominantly those involving the KK partners of singlet leptons) are typi- 
cally weaker than the ones for 7^7^ annihilation. Thus the contribution of coannihilation 
channels, despite that they do not suffer much from a Boltzmann suppression factor, 
Bi 0.9, is more than compensated by the increase in the effective number of degrees 
of freedom, see Eq. [13 This is in sharp contrast with supersymmetry where coannihi- 
lation generally decreases the relic density as the coannihilation channels typically have 
much larger cross section than the main channel. A large number of coannihilation pro- 
cesses contribute to the effective annihilation cross sections, primarily coannihilation with 
Higgs, '-f^H^ — )■ t, t{b) as well as a host of other self annihilation processes of the type 
IH^ XX, H^H^ XJC, l^H^ XJC, l^V^ -> XX', each individual channel con- 
tributing to a small fraction (less than 1%) of the total effective annihilation cross section. 
Here / stands for e,/i,r, stands for W^, and for any of the neutral or charged 
Higgs. Coannihilations processes involving quarks are completely negligible. Adding the 
loop-induced couplings between level 2 KK particles and SM particles has a significant 
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Figure 1: Vlh'^ as function of R'^ for m,, = 120 GeV, Ai? = 20 (left) and AR = 50 (right) 
including different processes as specified on the figure. Here 1-Ioop stands for one-loop 
couplings between level 2 and SM particles. The shaded region corresponds to the 3a 
preferred region obtained by WMAP [13]. 



impact on the relic density, see Fig. [H This is mainly because the new contribution from 
the process '-f^h^^ — > /i^"*" — )■ tb benefits from a resonance enhancement thus increasing 
significantly the effective annihilation cross section. This result depends very sensitively 
on the mass of the level-2 particle, a small downward shift in the mass, such as in the 
MUED model used in [11], where the renormalization scale is set to /i = 2R~^ for the 
level 2 masses, means that the pole effect is avoided at the LKP decoupling temperature. 
When including the contribution of h"^ and neglecting level 2 KK-particles in the final 
state, the prediction for the relic abundance is close to the one obtained including only 
annihilation processes. 

When allowing level-2 particles in the final state, mainly 7^ and h^,a'^,a^^, the relic 
abundance decreases sharply shifting the preferred value of the DM mass above the TeV 
scale. This is due to the important contribution of the coannihilation channels (/^7^ — )■ 
I'-f^) that are enhanced by the exchange near resonance of the n = 2 KK singlet lepton. 
Together these channels make up more than 50% of the (co) annihilation channels. As 
previously, other coannihilation channels each contribute to a small fraction of the total 
effective cross section. The contribution of the most important channels is illustrated in 
Fig. 2, where we have summed the contribution of all leptons in the initial states and all 
SM particles in the final state. Coannihilation channels involving lepton pairs contribute 
around 15% and their contribution is comparable to the one of Higgs channels 'f^H^ 
at large values of R~^. Contributions of the order of a few percent are found for the 
annihilation channels, 7^7^, as well as coannihilations of the type l^H^, H^H^ or 7^/^ into 
only SM particles. This still leaves around 10% contribution from all remaining channels, 
among these one finds notably channels involving gauge bosons such as V^H^ or V^l^ . 

The value of the cut-off scale A has an impact on the mass of the KK particles through 
logarithmic one-loop corrections, Eq. [TTl Increasing the scale to Ai? = 50 leads to heavier 
KK particles, in particular for KK lepton doublets and KK quarks, and has an impact on 
VLh"^. For example when ignoring the level 2 particles in the final state the contribution 
of coannihilation channels with KK leptons suffers from a larger Bolzmann suppression 
factor, this is partly compensate by an increase in the contribution of the pole (as 
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Figure 2: Relative contributions to the relic abundance at the LKP decoupling tempera- 
ture as a function of for rrih = 120 GeV, AR = 20. Here, summation over a class of 
initial states and all possible final states is performed with the exception of 7^/^ — )■ level 
2 SM which includes only processes with one level 2 particle in the final state and 7^/^ — > 
SM SM which includes only SM particles in the final state, stands for e^, ij,^,t^, u}, 
for W^, and for a^, h^^, h^. All remaining channels contribute less than 1%. 



well as the one of /i^, a^) and therefore in the contribution of the coannihilation channels 
involving Higgs. The net effect is an increase in Qh"^ by around 5%. On the other hand 
in the complete calculation including all channels, both the increase in the contribution 
of the /i^^ pole and the increase in the pole contribution of KK leptons which dominate 
the effective cross section lead to a 15% decrease of Qh"^. This implies an additional shift 
of close to lOOGeV in the preferred range for the DM mass. 

The value of the light Higgs mass which enters the loop corrections to the KK Higgs 
masses, has some impact on the prediction of the relic density. Increasing the mass from 
rrih = 120 GeV increases the contribution of coannihilation channels involving level 1 KK 
Higgses thus reducing see Fig. [31 The effect is noticeable only when nih > 200 GeV 
and is of the order of 7% for i?"^ = 1.3 TeV. For Higgs masses around nih ~ 220 GeV, 
the charged Higgs becomes the LKP [6]. 



4.1 Allowing for arbitrary mass splittings 

To illustrate the importance of the exact mass splittings between the KK states we gener- 
alize our model by allowing additional corrections to the fermion masses. In practice this 
is done by introducing a small shift in Zf, Eq. [HI This procedure guarantees that gauge 
invariance is maintained at the tree-level. This also implies that the relative mass shift is 
applied to both n = 1 and n = 2 KK leptons. First we illustrate the effect by modifying 
only the mass splitting of the KK partners of the lepton singlets, we assume generation 
universality. The free parameters of the model then include in addition to AR, R and rrih 
of the minimal model, five new parameters that affect the masses of the KK fermions, 

^Iri ^II^ ^QLi^dm Zu fi- 
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Figure 3: Contour plot of Qh"^ = 0.11 in the — rrth plane for KR = 20 (left) and 
AR = 50 (right). The shaded region correspond to the 3a WMAP range, 0.0952 < Qh'^ < 
0.1288, in the case where level 2 KK-particles in the final state are included (dark) or 
neglected (light grey). All coannihilation channels are taken into account. In the region 
above the full contour the LKP is the charged Higgs. 

We observe, see Fig.|4]for the case of the singlet leptons, that the value ofQh'^ increases 
with a smaller mass splitting. This might seem in contradiction with the discussion above 
since we had argued that the lepton coannihilation had the effect of decreasing Qh'^, as 
seen from Fig. 2, aO and cO. The main effect of a smaller mass splitting is to reduce the 
contribution of the channel o"(e}j.7^ — )■ 67^), indeed the resonance moves very near 
the threshold for the reaction and so does not contribute significantly to the thermally 
averaged cross section. This effect is more significant than the increase in the Boltzmann 
factor which can be at most 15% since in MUED for lepton singlets, B^j^ = 0.86 for 
R~^ = 1.3 TeV. In a sense the relic density is moving towards the value it would have 
if we had neglected the production of 7^ in the final state. Conversely an increase in 
the mass splitting leads to a mild decrease in fi/i^, here the Boltzmann suppression of the 
coannihilation channels is more than compensate by the decrease in the number of degrees 
of freedom. Note however that the relic abundance is insensitive to the mass splitting if 
it is more than 3%. 

We have also examined the effect of the mass splitting with the partners of the left- 
handed leptons. The effect follows the same trend although the influence on flh"^ occurs 
for splittings below 5%. The maximum increase in Qh^ is comparable to the one obtained 
for singlet leptons, see the left frame of Fig. |H Decreasing the mass of KK quarks on 
the other hand has the opposite effect as for leptons. A smaller mass splitting leads to 
a lower value for Qh"^, this is because in this case the factor Bi changes significantly and 
QCD processes of the type q^q^ qq give a large contribution. To illustrate this we 
consider the case where we shift the mass of the KK singlet d-type quarks, see the right 
frame of Fig. HI Finally we have also considered the implication of mass shifts for the 
Higgs. The KK-Higgs masses are modified either by increasing the light Higgs mass or by 
introducing a mass shift via the parameter Z^. In both cases this can lead to an increase 
of Qh'^ around 10% when the mass difference is a few per-mil. 

In summary keeping the mass splitting as a free parameter allows to find scenarios 
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that have Qh"^ in the range preferred by cosmological measurements for a lower KK scale. 
For example for = 1 TeV, we find Qh'^ = 0.11 when mgi^ — m^i = 1.7 GeV or 
TTigi^ — m^i = 7 GeV. On the other hand enlarging significantly the mass splittings between 
the LKP and all level 1 KK particles so as to reduce the contribution of the coannihilation 
channels would bring us back to the no coannihilation case with a preferred value for the 
mass of the DM candidate around SOOGeV, see Fig. [TJ 
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Figure 4: Qh"^ as function of a) Ag = (fn^i^ — m^i)/m^i (black) or Ag = {m^i^ — m^i)/m^i 
(red-dot) and b) A^^ = {rrij^i^ — m^i)/m^i for = 0.5, 1, 1.3 TeV and nih = 120 GeV. 
The other parameters are as in MUED. Blobs represent MUED points. The shaded region 
corresponds to the 3a preferred region obtained by WMAP. 



4.2 Dark matter searches 

The computation of the LKP-nucleon scattering cross section relevant for direct detection 
was performed with micrOMEGAs2.4 [27j. The elastic LKP-nucleon elastic scattering cross 
section depends both on the Higgs exchange as well as on the diagrams with exchange 
of level 1 KK-quarks. The amplitude is proportional to l/m1 for the first contribution 
and inversely proportionnal to the mass difference between the level 1 KK-quarks and the 
LKP, 1/A^ for the second contribution For the typical mass difference of the MUED 
scenario, A = m^i — m^i ~ 0.17m-yi the Higgs exchange is dominant. When computing 
the LKP-nucleon cross section we will use the micrOMEGAs option that includes the loop 
contribution to the KK quark exchange diagram. Although this contribution is computed 
exactly only for scalar particle exchange, it gives a better approximation than the tree- 
level calculation. Indeed at tree-level one can see the effect of the resonance contribution 
of the diagram when is such that rriq + m^i = m^i. This resonance effect is not 
physical and is just a sign that including the t-quark contribution by taking into account 
tree-level diagrams together with a coefficient for describing the t-quark content of the 
nucleon is not a good approximation [22] • To compute the elastic scattering cross section 
we use two sets of quark coefficients in the nucleon, the default values of micrOMEGAs , 
o'nN = 56, (To = 35 MeV, as well as values extracted from recent lattice calculations 
Ctj-m = 47, o"o = 42.9 MeV. The lattice calculations in particular indicate that the average 
value for the s-quark operator as = 50 MeV is lower than expected before [30] . 
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In MUED, cross sections are rather low, typically several orders of magnitude below 
the best limits of CDMS and Xenon. We compute the rescaled cross section on point-like 
nucleus [3T]. The rescaling factor takes into account the fact that 7^ does not account 
for all the dark matter. It is set to = f2/i^/0.0945 when DM is underabundant and to 
,^ = 1 otherwise. Although we display results for scattering on Ge^^, the protons and 
neutrons contribution do not differ much since the contribution from Higgs exchange is 
dominant, therefore aiKP-Ge ~ o'lkp-u ~ ctlkp-p- The rescaled cross section is rather 
stable around 1 x 10~^° pb, see Fig.[5]for the default parameters and is reduced by roughly 
a factor 2 using the lattice coefficients. 

In non-minimal models where the mass splittings are treated as free parameters, the 
LKP-nucleon scattering cross sections can increase by orders of magnitude due to the 
contribution of the quark exchange diagram which is proportionnal to 1/ (m^i — m^i)^ ^2E\ . 
To illustrate this we decrease the mass of the right-handed d quarks by treating as a 
free parameter. For example for = 1.3 TeV and Am ^ 7 GeV, we find a ~ 10^^ pb, 
about one order of magnitude below the limit of CDMS [32]. Of course, introducing such 
a small mass splitting has also an impact on the relic abundance, the new coannihilation 
channels with KK quarks reduce Qh"^ as discussed in section 14.11 Note that when the 
mass splitting becomes very small one can see the effect of the resonance contribution of 
the diagram when is such that rriq + m-yi = rriqi, this effect is present even using 
the loop improved calculation. 

For rrih = 220 GeV, ctlkp-n drops by almost one order of magntiude in the MUED 
case, this is because the Higgs contribution completely dominates. On the other hand 
the decrease in the cross section with the larger Higgs mass is more modest in the case 
where A^^ is small, see Fig. [51 This is because the Higgs contribution becomes sub- 
dominant. The spin dependent cross section is also typically small in MUED, for example 
al^P_p = 2.4 X 10-^pb and aff p_p = 2.4 x lO'^pb for R-^ = 1.3 GeV. These are almost 
six orders of magnitude below the current limits |331 EH ES] • 

For completeness we have also computed typical cross sections for KK-particles pro- 
duction at LHC with ^/s = 14 TeV . The largest cross sections are obtained for coloured 
particles. The dominant process are q^q^ pairs followed by q^g^, both with cross sections 
in the range 0{.l — 1) pb for the DM preferred mass scale, see Fig. O These results are 
in agreement with previous calculations [SH EH ESI EH] and are slightly suppressed as 
compared with the ones obtained using only tree-level KK masses. 

5 Conclusion 

We have performed a complete computation of the relic density of dark matter in the 
minimal UED model including all effects of level 2 particle in the intermediate and final 
state. We have shown that the production of 7^ in the final state reduces significantly the 
relic density thus shifting the preferred region from 500-600 GeV to well above the TeV 
scale (more precisely around 1.3TeV). The slight tension between the electroweak precision 
observables which favour R~^ > 600 GeV [2T1 [22] for a Higgs mass of 115 GeV and the 
dark matter observables is thus released. Indeed for R~^ = 1.3 TeV the contributions to 
the 5" and T parameters are suppressed, with in particular T ^ 10~'^. On the other hand 
the higher scale means that it will be harder to probe these scenarios at LHC as well as in 
direct detection. Furthermore because of the important contribution of the coannihilation 
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Figure 5: Rescaled LKP-nucleon cross section on Ge^^ vs rriLKP for "^h = 120 GeV, 
AR = 20 and 2 sets of quark coefficients ( {a^N, o'o) = (56 MeV, 35 MeV) (dash) or 
(47 MeV, 42.9 MeV) (dot) ) and for different values of the mass sphtting between the 
KK singlet d-quarks and the LKP including the MUED case (left panel). The MUED 
results for rrih = 220 GeV are also shown. In each line the region between the blobs is 
consistent with the 3a WMAP range. Rescaled LKP-nucleon cross section on Ge''^ vs rrih 



for R-^ = 1300 GeV, AR = 20 
with the 3a WMAP range. 



(right). In each line the region left of the blob is consistent 




R (TeV) 

Figure 6: a{pp — )■ q^q^,q^g^, g^g^,t^t^) vs R^^ at the LHC for AR = 20 using the 
CTEQ6M PDF's. The K-factor is not included. The shaded region corresponds to the 
3a preferred region obtained by WMAP. 



channels, the typical value of (av) relevant for indirect cross section signals is suppressed 
relative to the typical cross section expected for models that give Qh"^ = 0.1. Generalizing 
the model to allow for arbitrary mass shifts in the KK spectrum, we have shown that one 
could again increase the relic density so that agreement with WMAP was recovered for a 
LKP around the TeV scale in the case where the lepton NLKP were almost degenerate 
with the LKP. In this case the direct detection cross section could be strongly enhanced. 
We have also shown that these results not only depend sensitively on the mass differ- 
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ence between the level 1 particle and the LKP but also on the precise mass of the level 2 
particles that can be exchanged in s-channel in annihilation or coannihilation processes. 
In that sense making a precise theoretical prediction of the relic abundance of DM based 
on collider observables |10] in the event of the observation of KK-particles at the LHC is 
expected to be extremely challenging. Indeed it would require not only the measurement 
of the mass and couplings of the LKP but also a precise determination, in some cases 
better than the percent level, of the masses of level 2 particles. Such precision is not 
within the reach of the LHC especially for particles that are well above 1 TeV. 
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